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Abstract —One of the key steps in ultrasound image formation 
is digital beamforming of signals sampled by several transducer 
elements placed upon an array. High-resolution digital beam¬ 
forming introduces the demand for sampling rates significantly 
higher than the signals’ Nyquist rate, which greatly increases the 
volume of data that must be transmitted from the system’s front 
end. In 3D ultrasound imaging, 2D transducer arrays rather 
than ID arrays are used, and more scan-lines are needed. This 
implies that the amount of sampled data is vastly increased with 
respect to 2D imaging. In this work we show that a considerable 
reduction in data rate can be achieved by applying the ideas 
of Xampling and frequency domain beamforming, leading to a 
sub-Nyquist sampling rate, which uses only a portion of the 
bandwidth of the ultrasound signals to reconstruct the image. 
We extend previous work on frequency domain beamforming for 
2D ultrasound imaging to accommodate the geometry imposed 
by volumetric scanning and a 2D grid of transducer elements. 
We demonstrate high image quality from low-rate samples by 
simulation of a phantom image comprised of several small 
refiectors. We also apply our technique on raw data of a heart 
ventricle phantom obtained by a commercial 3D ultrasound 
system. We show that by performing 3D beamforming in the 
frequency domain, sub-Nyquist sampling and low processing rate 
are achievable, while maintaining adequate image quality. 

Index Terms —Array Processing, Beamforming, Compressed 
Sensing, Ultrasound 

1. Introduction 

Sonography is one the most widely used imaging modalities 
due to its relative simplicity and radiation free operation. 
It uses multiple transducer elements for tissue visualization 
by radiating it with acoustic energy. The image is typically 
comprised of multiple scanlines, obtained by sequential in- 
sonification of the medium using focused acoustic beams. 
Reflected signals detected at each transducer element are 
sampled prior to digital processing. Beamforming is the key 
step in image formation allowing for generation of receive 
sensitivity profile focused at any desired point within the 
image (2D) or volume (3D). The resulting beamformed signal, 
characterized by enhanced signal-to-noise ratio (SNR) and 
improved angular localization, forms a line in the image, 
which we refer to as beam. 

A. Motivation 

The aforementioned approach, used by most commercial 
systems today, is characterized by two important parameters 
- sampling and processing rate and frame- or volume-rate. 
Sampling rates required to perform high resolution digital 


beamforming are significantly higher than the Nyquist rate of 
the signal ID. Taking into account the number of transducer 
elements and the number of lines in an image, the amount of 
sampled data that needs to be transferred to the processing 
unit and digitally processed is enormous, even in 2D imaging 
setups, motivating methods to reduce sampling rates. In addi¬ 
tion, regardless of computational power, the frame/volume-rate 
in this approach is limited by the time required to transmit a 
beam, receive and process the resulting echoes, and to repeat 
the process for all image lines. 

Among the main focuses in the study of ultrasonic scanning 
is the development of real-time 3D ultrasound imaging, which 
overcomes major constraints of 2D imaging. 3D volume acqui¬ 
sition eliminates operator dependence in the imaging process 
- once the 3D data set is obtained, any plane within it is 
available for visualization by appropriate cropping and slicing. 
In addition, a variety of parameters can be measured from a 
3D image in a more accurate and reproducible way compared 
to 2D imaging O, O, (H. Moreover, many anatomical struc¬ 
tures, e.g. the mitral valve, are intrinsically 3D O, implying 
that their complex anatomy cannot be captured efficiently with 
2D techniques. 

A straight-forward approach to 3D volume acquisition is 
using a mechanically rotating ID probe IH. However, this 
technique suffers from extremely low volume rates, leading 
to unacceptable motion artifacts in echocardiography appli¬ 
cations. Fully sampled 2D arrays, an extension of the ID 
array to both lateral and elevation directions, are the most 
advanced technology for intrinsic 3D acquisition. Such arrays 
allow for significant improvement in frame rate and real-time 
capabilities. This is obtained by so-called ’parallel processing’, 
namely, electronically receiving data from several points in 
both lateral and elevation dimensions within the 3D volume 
simultaneously El, m. 

Being an optimal solution in terms of frame-rate, angular 
resolution and SNR, fully sampled 2D arrays pose several en¬ 
gineering challenges O, (91. Due to the significantly increased 
number of elements, which can be as high as several thousand, 
the main challenge from a hardware perspective is connecting 
the elements to electronic channels. In addition, the amounts 
of sampled data, acquired at each transmission, create a bot¬ 
tleneck at data transfer step and pose a severe computational 
burden on the digital signal processing hardware. To avoid too 
large connecting cables leading to unacceptable probe size and 
weight and to keep the electronics reasonable in power and 
size, as well as to reduce the data rates, numerous techniques 
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for element number reduction have been proposed. 

A straight-forward approach, referred to as sparse aperture, 
is to use only a subset of the 2D grid of elements upon recep¬ 
tion and/or transmission. Several studies investigate strategies 
for optimal subset choices 0, Ha-Ill which limit the 
reduction in image quality due to energy loss and high grating 
and side-lobes. In ifTSll Savord and Solomon present a sub¬ 
array beamforming approach allowing for significant reduction 
in the number of channels by sub-optimal analog beamform¬ 
ing, also referred to as micro-beamforming. This method was 
lately implemented in leading commercial systems. Another 
promising method, synthetic aperture, was adopted from sonar 
processing and geological applications (Tbl, IITtII . This ap¬ 
proach exploits multiplexing to control a fully sampled 2D 
array with a small number of electronic channels. Although 
providing improved image quality, synthetic aperture suffers 
from reduced frame rate and huge amounts of sampled data. 

Even when reducing the number of elements the amount 
of sampled data is still very large due to the high number 
of scan-lines. Consider ultrasonic imaging of a 3D volume, 
using K scan-lines in each one of K 2D cross-sections of the 
volume. Scanning the entire volume yields a total of K x K 
scan-lines. To maintain the angular resolution in each one of 
the K cross-sections in the 3D frame in comparison to 2D 
ultrasound imaging, one is forced to essentially quadrate the 
amount of data with respect to 2D imaging, given the same 
amount of transducers. 

B. Related work and Contributions 

In this work we present an approach for data rate reduction 
which can be applied in conjunction with any of the existing 
methods for element reduction. Our technique generalizes 
beamforming in frequency developed for 2D imaging. 

To achieve sampling and processing rate reduction in 2D 
Chernyakova and Eldar llT^ extended the concept of com¬ 
pressed beamforming IT9l and proposed performing beam¬ 
forming in the frequency domain. In this approach the Eourier 
coefficients of the beam are computed as a linear combination 
of those of the individual detected signals, obtained from 
their low-rate samples. When all the beam’s Eourier coeffi¬ 
cients within its bandwidth are computed, the sampling and 
processing rates are equal to the effective Nyquist rate. The 
beam in time is then obtained simply by an inverse Eourier 
transform. This approach is valid without any assumptions on 
the ultrasound signal structure. When further rate reduction 
is required, only a subset of the beam’s Eourier coefficients 
is obtained, which is equivalent to sub-Nyquist sampling and 
processing. Recovery then relies on an appropriate model of 
the beam, which compensates for the lack of frequency data. 

The work in mi demonstrates low-rate 2D ultrasound 
imaging, including the sub-Nyquist data acquisition step, low- 
rate processing and beamformed signal reconstruction. Low- 
rate data acquisition is based on the ideas of Xampling 
Ga-Gi, which obtains the Eourier coefficients of individual 
detected signals from their low-rate samples. More specifically, 
using Xampling we can obtain an arbitrary and possibly 
nonconsecutive set n, comprised of K frequency components. 


from K point-wise samples of the signal filtered with an 
analog kernel s*(t), designed according to n. In ultrasound 
imaging with modulated Gaussian pulses the transmitted signal 
has one main band of energy. As a result the analog filter takes 
on the form of a band-pass filter, leading to a simple low-rate 
sampling scheme mi. The choice of n dictates the bandwidth 
of the filter and the resulting sampling rate. 

In 3D imaging the same low-rate sampling scheme can be 
applied to the individual signals detected at the elements of the 
2D transducer, leading to considerable rate reduction, as elab¬ 
orated on in Section IIII-Bl However, to benefit from the rate 
reduction, 3D beamforming must be performed in frequency 
similarly to the 2D setup. We prove that the relationship 
between the beam and the detected signals in the frequency 
domain, the core of beamforming in frequency, holds in the 
3D imaging setup as well 

In this work we derive a frequency domain formulation 
of beamforming that accounts for the 2D geometry of the 
transducer array and the 3D geometry of the medium. We 
show that, similarly to 2D imaging, 3D frequency domain 
beamforming (EDBE) can be implemented efficiently due to 
the decay property of the distortion function translating the 
dynamic beamforming time delays into the frequency domain. 
When sub-Nyquist sampling and processing are applied, signal 
structure needs to be exploited to recover the beam from the 
sub-Nyquist set of its Eourier coefficients. To this end we prove 
that a 3D beamformed signal obeys a finite rate of innovation 
(ERI) 1^ model, just as in 2D. 

We next report the results of simulations and experiments 
verifying the performance of the proposed method in terms of 
the lateral point spread functions (LPSE), axial point spread 
functions (APSE) and SNR. Einally we incorporate 3D EDBE 
to a commercial imaging system performing analog sub¬ 
array beamforming and show that the above techniques are 
compatible, namely, 3D EDBE does not introduce additional 
image degradation. 

The rest of the paper is organized as follows: in Section JIl 
we review standard time-domain processing for a 3D imaging 
setup. In Sections Hill and HVl we describe the principles of 3D 
EDBE, image reconstruction and the achieved rate reduction. 
In Section [V| the results and comparison to time-domain 
beamforming are presented. Section [yl| concludes the paper. 

II. Implementation of Beamforming in Time 

Beamforming, a basic step required by all ultrasound based 
imaging modalities, is a common signal-processing technique 
that enables spatial selectivity of signal transmission or recep¬ 
tion ED- In ultrasound imaging it allows for SNR and lateral 
resolution improvement. Modern imaging systems transmit 
and receive acoustic pulses using multiple transducer elements. 
These elements comprise an array, generating a transmitted 
beam which is steered spatially by applying appropriate time 
delays to each element. The transducer receives acoustic pulses 
scattered by tissue structures, which are then sampled and pro¬ 
cessed digitally to reconstruct an image line. Reconstruction is 
performed with a technique known as dynamic beamforming, 
where the image quality is enhanced by summing the signals at 




3 


individual elements after their alignment by appropriate time- 
delays. 

To derive frequency-domain implementation of 3D beam¬ 
forming we begin by introducing standard time-domain pro¬ 
cessing. Consider a grid of M x transducers located in 
the x-y plane, depicted in Fig. [T] The geometry imposed by 
3D ultrasound imaging requires the use of two steering angles 
and thus a 2D array of transducers. The entire grid transmits 
pulses into the tissue. We note that the grid may have a small 
curvature along the 2 ; axis, so the array elements do not lie 
in the same plane. For the sake of simplicity, this type of 
curvature is not displayed in Fig. [T] 



z 

Fig. 1: M X transducers placed in the x-y plane. An acoustic 
pulse is transmitted in a direction Ox^Oy. The echoes scattered 
from perturbations in the radiated tissue are received by the 
array elements. 


We choose a reference element, (mo, no), placed at the 
origin, and denote the distances along the x and y axes to 
the {m^n) element by 6m respectively. We also denote 
the height of the (m, n) element with respect to the origin by 
6^ ^ for the case where there exists a curvature along the ^ 
axis. Note that we assume 6^^ =0, so that the reference 

element is not necessarily included in the transducer grid, and 
is defined for mathematical convenience. Let us consider a 
pulse transmitted along a scan-line specified by spatial angles 
Ox^Oy. Setting t = 0 at the moment of transmission from the 
(mo, no) element, it can be shown that at time t > 0 the pulse 
reaches the coordinates: 


with 


z{t)) = ct{xe,ye, ze), 
sin Or cos O.j 


xe = 


ye = 


Ze = 


1 — sin^ Ox sin^ Oy 
cos Ox sin 0.J 


0., 


1 — sin^ Ox 
cos Ox cos Oy 


1 — sin^ Ox sin^ Oy 


( 1 ) 


( 2 ) 


Here c is the propagation velocity in the medium. A point 
refiector located at this position scatters the energy, such that 


the echo is detected by all array elements at a time depending 
on their locations. 

Denote by (^rn,n(^; Oy) the signal detected by the (m, n) 
element and by rm,n{t] Ox^Oy) the time of detection. Then: 

= (3) 

where 


dm,n{t'’) Ox-! Oy) — (4) 

^{x{t) - (5„)2 + {y{t) - (5„)2 + {z{t) - 

is the distance traveled by the refiection. Beamforming in¬ 
volves summing the signals detected by multiple receivers 
while compensating for the differences in detection time. 
Using ©, the detection time at (mo, no) is 

Ox^ Oy) — 2f since 6mo — 6fiQ — ^m,n — 

We wish to apply a delay to (pm,n{t]0x^0y) such that the 
resulting signal, denoted by 0m,n{t;Ox^Oy), satisfies: 


0xi Oy) — ^m,n(nm,n(^5 Oxi Oy)] Ox-, Oy). 


Doing so, we can align the refiection detected by the 
{m^n) receiver with the one detected at (mo, no). Denoting 
Tm,n{t;0x,9y) = fm,n{t/2-,9x,0y) and using the follow- 
ing aligned signal is obtained: 

Ox', Oy) — rn,nij~rn •) Oxt Oy)) Ox) Oy)^ 

Ox) Oy) = ( 5 ) 

^ 4|7m,nP - {'JmXe + 'Jnye + lm,nZe)) , 

where we defin ed 7 ^ = 6m/c, 7 n = <5„/c, „ 

and \jm,n\ = + (7m,n)^- 

The beamformed signal may now be derived by averaging 
the aligned signals. We assume that the echo reception process 
involves a subset of the transducer array, denoted by C 

{(m, n)| 1 < m < M, 1 < n < A^}: 

^{t']0x)0y) = — \ ^ Ox) Oy) . ( 6 ) 

iVRX — 

(m,n)eM 

Here A^rx=|A^| is the number of transducers participating 
in the reception process. We note that in order to obtain 
optimal performance in terms of SNR and angular resolution, 
all transducer elements should be used. However, as mentioned 
in Section U the number of active elements is often reduced 
due to practical constraints. 

The beamforming process is carried out digitally, rather than 
by manipulation of the analog signals. The signals detected at 
each element must be sampled at a sufficiently high rate to 
apply high-resolution time shifts defined in (O. This implies 
that the signal is sampled at rates significantly higher than its 
Nyquist rate, in order to improve the system’s beamforming 
resolution and to avoid artifacts caused by digital implementa¬ 
tion of beamforming in time. From now on we will denote this 
rate as the beamforming rate fs, which usually varies from 4 
to 10 times the transducer central frequency [HI, 

We conclude this section by evaluating the number of 
samples typically required to obtain a single volume for 
some predefined image depth. Our evaluation is based on the 
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imaging setup used in the simulation displayed in Section 
The simulation assumes an ultrasonic scanner comprising a 
32x32 grid of transducers, all of which are active both on 
transmission and reception (A^rx = 1024). Such an array 
constitutes a reference for comparison of image quality result¬ 
ing from different methods for data rate reduction. The radial 
depth of the scan is set as r = 5.5 cm with a speed of sound 
of c = 1540m/sec, yielding a time of flight of T = 2r/c 
71.43 /isec. The acquired signal is characterized by a band¬ 
pass bandwidth of 1.4 MHz centered at a carrier frequency of 
/o = 3 MHz. It is sampled at a rate of fs = 18.25 MHz 
to provide a sufficient beamforming resolution leading to 
N = 1304 samples taken at each transducer. Every frame 
contains 21 x 21 scan-lines, such that the scanned volume is 
a square pyramid with an opening angle of 14.3°. This set of 
scanning angles is a relatively narrow set with a typical margin 
between subsequent beam lines. Therefore, assuming that it is 
possible to sample all 1024 elements to obtain optimal image 
quality, the total number of samples that must be processed to 
display a single frame is 21 x 21 x 1024 x 1304 = 5.89 • 10^. 
This number of samples is huge even for a moderate imaging 
depth of 5.5 cm; the imaging depth typically required for 
cardiac imaging is around 16 cm. Achieving a reasonable 
frame-rate using such an amount of samples is infeasible for 
any low-cost ultrasound machine. Therefore, even assuming a 
hardware solution allowing for connection of all the transducer 
elements to electronic channels, the amount of data is still a 
bottleneck. As a result sparsely populated arrays of transducer 
elements may be used. This typically causes a reduction 
in angular resolution and, more significantly, low SNR. A 
solution that reduces the amount of samples while using the 
entire transducer grid in the reception stage will address this 
problem. 


HI. Beamforming in Frequency 

To substantially reduce the number of samples taken at 
each transducer element we aim to use the low-rate sampling 
scheme proposed in ifT^ . To this end we derive a frequency- 
domain formulation of 3D beamforming allowing to compute 
the Fourier coefficients of the beam from the detected signals’ 
low-rate samples. In this section we show that similarly to 
2D imaging the Fourier coefficients of the 3D beam can be 
computed as a linear combination of the Fourier confucianists 
of the received signals. We note that due to the dynamic nature 
of beamforming, such a relationship is not trivial and requires 
appropriate justification. 

A. Beamforming in Frequency Scheme 

We start from the computation of the Fourier series co¬ 
efficients of the beamformed signal ^{t]0x^0y). It is shown 
in Appendix [B] that the support of ^{t^Ox^Oy) is limited to 
[0 ,Tb( 6>^,6>^)), where TsiOx^Oy) is given by: 

TB{0x,0y)= mm T-l^{T;0x,Oy), (7) 

{m,n)EA 4 

with rm,n{F,0x^0y) defined in (O. It is also shown that 
TsiOx^Oy) < T, where T is defined by the transmitted pulse 
penetration depth. 


Consider the Fourier series of the beamformed signal, 
{c[k]}k, in the interval [0,T]: 

= ^{t\ex,ey)I^QXB(e.,ey))e~"^'"^dt, (8) 

where fa^h] is the indicator function, plugged in to cancel 
noise since the useful information in ^{t]0x^0y) is restricted 
to f)^TB{0x^0y)).ln order to find a relation between c[k] and 
the Fourier coefficients of we substitute (|6]) 

into dSl): 


2 _^ 2 

(m,n)GA4 ^ 


Nrx . . 

where we defined 


(9) 


1 

-L Jo 


Substituting the integration variable u with 

r = rm,n{u] Ox, Oy) we get 


- h 1^ 

^ _ 'I /m,n I 

r - {imXe + 'ynVe + lm,nZe) ’ 

, 1-2 + - 2t • {'ymXe + 7n2/0 + lm,n^e) , 

du = - 2 "^-ar, 

[r - {-fmXe + InVe + 7m,n^e)] 

where xo^ye, zq are defined in Q. Plugging this into ([TOb and 
renaming the integration variable r ^ t, result in 


1 

^Ak] - y 


Ox, Oy^^ri 


with 


i{t;0x,0y)e 

( 11 ) 


(lk,m,n{t]0x,0y) — I[\j^^^\,Tm,n(TBi9x,9y)-,9x,ey)){t)x 

+ |7m,nP - • hmXe + + 7m 

- ^^ -^x 

(t - {yrnXe + 7n^0 + 7m,n^^)) 

f a • (7m^0 +7n^0 +7m,n^0) - l7m,np 

exp < —i — k - j ^-r- 

[ ^ V ^ + lny9 + 7m,n^0j 


( 12 ) 



Note that in contrast to (Uni), (dUi contains a non-delayed 
version of ipm,n{t;0x,0y), while the delays are applied 
through the distortion function qk,m,n{t',0x,0y), defined in 
dEl). This allows us to express (fm,n{t’,0x,0y) in terms of 
its Fourier series coefficients, denoted by Cm,n[l]- We also 
make use of the Fourier coefficients of qk,m,n{F,0x,0y) with 
respect to [0,T], denoted by Qk,m,n]9x,9y[l], and rewrite ([T71) 
as follows: 


1 

/ Qk,m,nA0xAy)e~^^^’"~’'At 

I ^ -^0 

~ ^ ^ ^]Qk,myn]0x,Oy [^]* 


( 13 ) 
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The substitution of the distortion function by its Fourier coeffi¬ 
cients effectively transfers the beamforming delays defined in 
(O to the frequency domain. We note that qk,m,n{t] Ox^Oy) is 
independent of the received signals, namely, it is defined solely 
by the array geometry. Its Fourier coefficients, therefore, are 
computed off-line and stored as a look-up-table (LUT). 

According to Proposition 1 in C3, which can be easily ex¬ 
tended to the 3D imaging setup, Cm,n[k] can be approximated 
sufficiently well when we replace the infinite summation in 
(fTSl) by a finite one: 

L2 

Cm,n ,m,n-,9x,0y [/]. (14) 

l=-Li 

The Fourier coefficients of the beam, c[k], can now be easily 
calculated by plugging ([T4l) into (O: 

1 

E E ^m,n [k - l]Qk ,m,n-,9x-,0y [«]• (15) 

(m,n)EAi l= — Li 

The approximation in (O relies on the decay properties of 
{Qk,m,n]9x,9ym}' According to the results reported in ifT^ 
most of the energy of the Fourier coefficients of the 2D 
distortion function is concentrated around the DC compo¬ 
nent, allowing for efficient implementation of beamforming 
in frequency. This decaying property is retained in 3D beam¬ 
forming: numerical studies show that most of the energy of 
{Qk,m,n-,9x,9ym} Concentrated around the DC component, 
irrespective of the choice of k^m^n^Ox^Oy. We assume that 
for I < -Li and I > L 2 , {Qk,m,n-,9x,9y[l]} are several orders 
of magnitude lower and thus can be neglected. The choice of 
L 1,1/2 controls the approximation quality. We display these 
decay properties in Fig. [2l where [^] is plotted as 

a function of I for k = 300, m = 7, n = 21, Ox = 0.28 [rad] 
and Oy = 0.36 [rad]. 



Fig. 2: {Qk,m,n-,9x,9y[l]}, the Fourier coefficients of 

qk,m,n{t]0x^0y), display a rapid decay rapidly around 
the DC component. 

B. Rate Reduction by Beamforming in Frequency 

We now show that FDBF allows one to generate a frame 
using a reduced number of samples of the individual signals 


with respect to time-domain beamforming. When the signal’s 
structure is not considered, this is done by avoiding the 
oversampling factor required by digital implementation of 
time-domain beamforming. In this case the processing is 
performed at the effective Nyquist rate defined by the signal’s 
effective bandwidth. Further rate reduction can be obtained 
when the FRI structure of the beamformed signal is taken into 
account and compressed sensing (CS) techniques are used for 
its recovery (231, (221 • 

As can be seen in (Hsi), in order to calculate an arbitrary 
set n of size /C of Fourier coefficients of the beamformed 
signal, only JC F Li + L 2 Fourier coefficients of each one of 
the individual signals are required. The image line is then re¬ 
constructed from the beamformed signal’s Fourier coefficients 
{c[k]}. Calculating the entire set of Fourier coefficients in the 
bandwidth of the beamformed signal (3, |/3| = B, implies 
B » I/i + 1/2 and, therefore, allows one to obtain all the 
information in the frequency domain while avoiding oversam¬ 
pling required by time-domain beamforming. This is due to the 
fact that the low-rate sampling scheme described in Section iFB] 
obtains BFL 1 FL 2 ~ B Fourier coefficients of the individual 
signals required for FDBF from their B low-rate samples. 
Thus, performing FDBF by calculating the entire bandwidth of 
the beamformed signal achieves an approximately N/B rate 
reduction factor with respect to time-domain beamforming, 
where B is the number of Fourier coefficients in the bandwidth 
of the beamformed signal and N is the number of samples 
required by the beamforming rate fs. 

Further rate reduction is possible by acquiring a part of 
the bandwidth of the beamformed signal, /i C /3, |/i| = M. 
We may calculate it from M + I /1 + I /2 ~ M samples of 
the individual signals, which are sampled at a rate that is 
N/M lower than the standard beamforming rate fs. In Section 
nYi we take advantage of the beamformed signal structure 
to reconstruct the beam from its partial frequency data. A 
detailed discussion on the achieved rate reduction is provided 
in Section IV^ 

IV. Recovery Method from Sub-Nyquist Samples 

When all the beam’s Fourier coefficients within its effective 
bandwidth are computed, the beam in time is recovered by an 
inverse Fourier transform. When only a subset of the coef¬ 
ficients is obtained by sub-Nyquist sampling and processing, 
we exploit the structure of the beam to reconstruct it from its 
partial frequency data. 

According to (20|, we can model the detected signals at the 
individual transducer elements, {(Pm,n{t]0x,0y)}(^rn,n)eM^ 
FRI signals. That is, we assume that the individual signals 
can be regarded as a sum of pulses, all replicas of a known 
transmitted pulse shape: 

L 

^y) — ^ ^ ( 16 ) 

1=1 

Here h{t) is the transmitted pulse shape, L is the number 
of scattering elements in the direction of the transmitted 
pulse {Ox^ Oy), are the unknown amplitudes of the 

refiections and {ti^rn,n}f=i are the times at which the refiection 
from the l-th scatterer arrives at the (m, n) element. 
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It is shown in Appendix that the beamformed signal in 
3D imaging approximately satisfies the FRI model, just as it 
does in 2D imaging ca. Namely, it can be written as 

L 

(17) 

1 = 1 

where h{t) and L are defined as above, are the 

unknown amplitudes of the refiections and are the 

times at which the reflection from the l-th scatterer arrives 
at the reference element (mo, no). 

Having acquired the Fourier coefficients c[k] as described 
in the previous section, we now wish to reconstruct the 
beamformed signal. Since the beam satisfies the FRI model 
our task is to extract the unknown parameters, and 

{U}iLi that completely describe it. 

Using ([TtI) the Fourier coefficients of Ox^ Oy) are given 
by: 



L 

= h[k] ^ , (18) 

1=1 

where h[k] is the k-th Fourier coefficient of h{t). 

By quantizing the delays with quantization step 

Ts = such that ti = qiTg for qi G Z, we may write 
the Fourier coefficients of the beamformed signal as: 

AT-l 

c[k\ = h[k\ (19) 

1=0 

with N = [T/TgJ, bi = bi^i^qi and 6a,b is the Kronecker delta. 
We conclude that recovering the beamformed signal in time 
is equivalent to determining bi in ([T9l) for 0 < / < — 1. In 

vector-matrix notation ([18]) cab be rewritten as: 

c = HDb = Ab, (20) 

where c is a vector of length 1C with k-ih entry c[k], H is a 
1C ^ 1C diagonal matrix with k-i\\ entry h[k], D is a. 1C x N 
matrix whose rows are taken from the A x A DFT matrix 
corresponding to the relevant Fourier indices of ^{t;0x,0y), 
and b is a column vector of length N with l-th entry bi. 

We wish to extract the values of b, which fully describe the 
beamformed signal. To do so, we rely on the assumption that 
a typical ultrasound image is comprised of a relatively small 
number of strong reflectors in the scanned tissue. In other 
words, we assume the vector b to be compressible, similarly 
to lIT^ . We then find b by solving an optimization problem: 

min||b||i s.t. ||Ab - c ||2 < (21) 

b 

In practice, we solve (121]) using the NESTA algorithm ll^ 
which works well when the signal of interest has high dy¬ 
namic range. NESTA uses a single smoothing parameter, /i. 


selected based on a trade-off between accuracy and speed of 
convergence. We choose this parameter empirically to achieve 
optimal performance with respect to image quality. 

To summarize this section, a step-by-step description of the 
3D low-rate imaging process is given in Algorithm [T] 


Algorithm 1 Image acquisition algorithm 

1: Calculate the Fourier coefficients of qj,m,n {t]0x^0y), de¬ 
fined in (Ha. This calculation is performed offline and 
does not affect the system’s real-time performance. 

2: Choose the approximation quality by determining 
defined according to the decay properties of 
{Qk,m,n;e^,ey [(]};, displayed in Fig. H An adequate ap- 
proximation can be performed by choosing I/i,I /2 to be 
no greater than 10. 

3: Choose a subset hz of Fourier coefficients of the beam- 
formed signal to be used in reconstruction. 

4: Acquire the Fourier coefficients of the individual signals 
relevant for reconstruction from low-rate samples, accord¬ 
ing to m. At each transducer element {m,n) G M, 
{cm,n where /ci, k 2 are the lowest and highest 

indices in the subset. 

5: Perform the calculation in ([14]) . 

6: Compute the beamformed signal’s Fourier series coeffi¬ 
cients: 

~ bm,n[k]. 

im^M 

7: Solve the optimization problem (121]) to extract the vector 
b that characterizes the beamformed signal. 

8: Incorporate the known temporal shape of the pulses, h (t), 
onto the vector b, and perform standard postprocessing 
steps, such as log-compression and interpolation. 


V. Simulations and Results 

To analyze the performance of the outlined methodology 
relative to standard time-domain beamforming in a manner 
independent of the specifics of any individual system, a k- 
Wave IZTll simulation of a 3D ultrasound system is presented. 
We first simulate the acoustic imaging of a noise free volume 
containing three point scatterers and analyze the effect of 
the achieved rate reduction on lateral and axial point spread 
functions. The performance of low-rate 3D FDBF is compared 
to that of standard time-domain processing. As mentioned in 
Section H] the number of samples can be reduced when only a 
partial set of transducer elements is used upon reception. To 
compare this approach to the proposed method, time-domain 
beamforming is also performed using the data collected only 
from the elements placed along the array’s main diagonals. 
We next show the advantage of 3D FDBF in the presence of 
noise compared to time-domain beamforming with the reduced 
number of elements. Finally we show that the proposed ap¬ 
proach can be incorporated into a commercial imaging system 
performing sub-array beamforming to reduce the number of 
channels. The results verify that rate reduction obtained by 3D 
FDBF does not introduce further image quality degradation. 
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A. Simulation Setup 

We simulate acoustic imaging of a volume of size 28 mm x 
28 mm x 55 mm. The volume contains three point reflectors, 
placed at depths 26, 31.5 and 37 mm from the center of a 
square planar 2D transducer grid. The reflectors are located 
around Ox = —7.5°, 0° and 7.5°, respectively, with Oy = 
0°. The reflector at depth 31.5 mm is located at the focus 
point of the transmitted pulse. The array is comprised of 32 x 
32 = 1024 elements, spaced 140 pm apart. The central pulse 
frequency is 3 MHz with a bandwidth of 1.4 MHz and the 
sampling rate is fs = 18.25 MHz. A penetration depth of 
T = 2r/c 71.43 /isec yields N = 1304 samples at each 
transducer element, so that the bandwidth of the beamformed 
signal contains B = 200 Fourier coefficients. A single volume 
comprises 21 x 21 scanned angles. 

Denote by n of cardinality K the set of Fourier coefficients 
of the beam, obtained by the proposed method. To verify the 
performance for different rate reduction factors the collected 
data is processed using our technique with 1C = B, 1C = B /2 
and 1C = Bl‘i corresponding to the entire, half and one third of 
the effective bandwidth respectively. The results are compared 
to those obtained by time-domain beamforming performed 
both using full and diagonal grids upon reception. 

First, we assess the amount of samples required to obtain a 
volume. For time-domain beamforming using the full grid, we 
must process 21 x 21 x 1024 x = 5.89 • 10^ samples. Using 
only the main diagonals of the transducer grid, this amount 
reduces to 21 x 21 x 64 x = 36.81 • 10^ samples. Applying 
FDBF, the reconstruction relied on /C = 200, 1C = 100 and 
/C = 67 Fourier coefficients of the beamformed signal. To 
calculate these coefficients, as described in Section luD we 
chose Li = L 2 = 10. The total amount of Fourier coefficients 
required at each transducer is z/ = /C-f-Li +L 2 . A mechanism 
proposed in |[T^ allows us to obtain these coefficients from u 
samples of individual detected signals. Thus, a single volume 
can be produced by processing a total of 21 x 21 x 1024 x 
u samples. The total number of samples required for each 
processing method is displayed in Table U We note that FDBF, 
using about half the bandwidth of the beam, is comparable to 
time-domain processing using only the diagonal elements of 
the grid in terms of processing rate. 

Cross-sections of the resulting 3D volumes are displayed in 
Fig. [ 3 ] It can be seen that all reflectors are clearly seen for 
all processing methods displayed. We note that the frequency- 
domain beamformed images display lower noise levels than 
the time-domain beamformed image reconstructed using a 
partial set of the transducer grid. The advantage of FDBF in 
the presence of noise is discussed in detail in Section IV-CI 

B. Lateral and Axial Point Spread Functions 

We next compare our proposed method to time-domain 
beamforming by calculating the LPSF and APSF character¬ 
izing each processing method. The LPSFs are acquired for 
the reflector placed at the transmit focus point and plotted on 
constant-r arcs. The APSFs show the sum of the constant-^a^ 
and Oy lines on which point reflectors are placed. The LPSFs 
are normalized such that the maximum at Oy = 0° is set to 1, 


Processing method 

Number of samples 

Time 

Full grid 

5.89 • 10® 

Time 

Diagonal grid 

36.81 • 10® 

Frequency 

B coefficients 

99.8 • 10® 

Frequency 

B/2 coefficients 

54.64 • 10® 

Frequency 

B/3 coefficients 

39.74 • 10® 


TABLE I: Number of samples per volume for each processing 
method. 

while the APSFs are normalized to unit energy. The PSFs are 
presented in Figs. [4] and [5l 


- Time-domain beamforming - full grid 

- Time-domain beamforming - diagonal grid 

-Frequency-domain beamforming - B DFT coefficients 

- Frequency-domain beamforming - Bj^ DFT coefficients 

- Frequency-domain beamforming - Bj?* DFT coefficients 


Normalized Lateral PSF 



^y[°] 


Fig. 4: Normalized LPSFs for various processing methods. 


Processing 

method 

Full width at 
half maximum 

Average 

side-lobes 

First side- 
lobe’s peak 

Time 

Full grid 

1.50° 

-29.63 dB 

-22.93 dB 

Time 

Diagonal grid 

1.51° 

-28.53 dB 

-21.40 dB 

Frequency 

B coefficients 

1.51° 

-29.01 dB 

-22.70 dB 

Frequency 

B/2 coefficients 

1.52° 

-29.37 dB 

-22.90 dB 

Frequency 

B/3 coefficients 

1.47° 

-29.47 dB 

-22.74 dB 


TABLE II: LPSFs properties. 


The properties of the LPSFs are displayed in Table HIl We 
see that the LPSFs obtained with 3D FDBF for different 
rate reduction factors display very similar properties to the 
LPSF acquired with time-domain beamforming using the 
full transducer grid, and exhibit improved results over the 
LPSF acquired with time-domain beamforming using only the 
diagonals of the grid. This is an expected result, since our 
method reconstructs the axial lines of the image and does not 
have a direct effect on the lateral resolution. The widths of the 

























(a) 


(b) 


(c) 


(d) 


(e) 



(f) (g) (h) (i) G) 


Fig. 3: Cross-sections of the simulated 3D imaging of three point reflectors placed on a plane, (a)-(e) display the Oy = 0° 
cross-section, while (f)-(j) display the Ox = 0° cross-section, (a), (f) display images acquired with time-domain beamforming, 
using the entire transducer grid, (b), (g) display images acquired with time-domain beamforming, using the diagonals of the 
transducer grid, (c) and (h), (d) and (i), (e) and (j) display images acquired with FDBF using B.Bjl and 5/3 DFT coefficients 
of the beamformed signal, respectively. 


- Time-domain beamforming - full grid 

- Time-domain beamforming - diagonal grid 

-Frequency-domain beamforming - B DFT coefficients 

- Frequency-domain beamforming - B/2 DFT coefficients 

- Frequency-domain beamforming - Bj?) DFT coefficients 



2 2.5 3 3.5 4 4.5 5 

2 ; [cm] 


Fig. 5: Normalized APSFs for various processing methods. 


peak located at the focus point, acquired for each reconstructed 
method, are shown in Table [Till As seen in Table [III| and 
Fig. E the APSFs deteriorate when the amount of Fourier 
coefficients used to reconstruct the image is decreased. We 


Processing method 

Full width at half maximum 

Time 

Full grid 

0.084 mm 

Time 

Diagonal grid 

0.082 mm 

Frequency 

B coefficients 

0.086 mm 

Frequency 

B/2 coefficients 

0.107 mm 

Frequency 

B/3 coefficients 

0.138 mm 


TABLE III: APSFs properties. 


note that energy leakage from the peaks is increased when 
less Fourier coefficients of the beamformed signal are used 
in the reconstruction process. However the effect of APSE 
deterioration becomes visible only when less then half the 
bandwidth is used for signal reconstruction. 

Reducing the amount of transducer elements enhances noise 
levels in the image and deteriorates the lateral resolution, 
while our proposed reconstruction method affects mainly the 
axial resolution. Acknowledging this fact, we may consider 
a midway approach, where rate reduction is achieved both 
by reducing the amount of transducer elements and applying 
frequency-domain beamforming. The more dominant factor 
for rate reduction will be dictated according to the trade-off 
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Stated above. In addition to the axial and lateral resolution, 
another important feature that has to be regarded is the SNR. 
It can be seen in Fig. [5] that the line acquired with time-domain 
beamforming using partial grid data contains high levels of 
noise at the far-held, since less transducer elements participate 
in the delay-and-sum process. This holds even when no noise 
was incorporated in the simulation - the “effective noise” stems 
from rehections corresponding to the numerical solution of the 
simulation code. 

C. Simulation with Noise 

To show the advantage of our proposed method over the 
partial grid time-domain beamforming method in terms of 
SNR, another simulation is conducted. A pulse is transmitted 
in the = 0, Oy = ^ direction. A single large rehector is 
placed at the focus depth of the transmitted pulse. The signals 
detected at the transducer elements are contaminated by white 
Gaussian noise imitating the thermal noise of the system. We 
then proceed to reconstruct the = 0, Oy = ^ beam using 
all hve methods described in Section IV-AI In addition, clean 
beams, without the addition of noise, are obtained for all hve 
methods. We denote the noisy and clean beamformed lines by 
^ noise (t) and ^ciean(^), respectively. 

We dehne the SNR of a beam as the ratio between the 
energy stored in a segment of length 5 A around the main peak 
of the beam, where A is the wavelength corresponding to the 
carrier frequency of the transmitted pulse, and the energy of 
the noise in the beamformed line, dehned as n{t) = ^noise(^) — 
^clean(^)- That iS! 

SNR 101ogio(^ J\n{t)\^dt )' 

The results are displayed in Table [Ivl As expected, the reduc¬ 
tion in number of elements participating leads to a dramatic 
reduction in SNR. In contrast, 3D FDBF displays higher SNR 
even over the time-domain processing when the entire grid of 
transducer elements is used. A remarkable point is that the 
SNR increases when less Fourier coefficients are involved in 
the reconstruction of the beam. This is not surprising since the 
noise is equally spaced over the entire spectrum of the system 
- the fewer Fourier coefficients used in the reconstruction 
process, the less noise involved. 

To conclude, 3D FDBF using half the bandwidth is compa¬ 
rable to standard time domain processing with diagonal trans¬ 
ducer elements in terms of data rate reduction, but outperforms 
it by 20 dB in terms of SNR. 


Processing method 

SNR 

Time 

Full grid 

20.83 dB 

Time 

Diagonal grid 

8.92 dB 

Frequency 

B coefficients 

26.69 dB 

Frequency 

B/2 coefficients 

28.68 dB 

Frequency 

B/3 coefficients 

29.48 dB 


TABLE IV: SNR of processing methods. 


D. Application on a Commercial System 

We now demonstrate our method on data collected using a 
commercial 3D ultrasound system, while imaging a phantom 
of a heart ventricle. Shown in Fig. [6] are images of the entire 
volumetric scan, taken from a specific angle, for demonstra¬ 
tion. The frames are reconstructed using time-domain beam¬ 
forming and frequency-domain beamforming with /C = 5/2 
coefficients of the beamformed signal. The complex structure 
of the phantom allows us to test the performance of the pro¬ 
posed method on volumes containing multiple strong reflectors 
as well as weak refiectors known as speckle. 

The transducer grid is comprised of 2000 transducers. The 
entire grid participates in the transmission stage, while analog 
sub-array beamforming is performed in the reception stage. 
This sub-optimal processing method is required to adjust the 
number of elements to the number of electronic channels. We 
processed the collected data in the same manner as in the 
previous section, using time-domain beamforming requiring 
3120 samples per image line, and frequency domain beam¬ 
forming for /C = 5 and 1C = B/2 with B = 506. When B 
Fourier coefficients are computed, the processing is performed 
at the signal’s effective Nyquist rate and oversampling is 
avoided leading to 6 fold rate reduction. FDBF with B/2 
coefficients implies 12 fold rate reduction leading to sampling 
and processing at a sub-Nyquist rate. Cross-sections of the 3D 
frames acquired are displayed in Fig. [71 

The image obtained by low-rate FDBF for /C = 5 is 
virtually identical to one obtained by standard time-domain 
processing at a high rate. This result is expected since for 
1C = B Sill the information is obtained in frequency. We 
also note prominent similarity between the image obtained 
at the sub-Nyquist rate and the original one. In particular, the 
strong refiectors and the speckle pattern are preserved. The 
above results prove that FDBF can be combined with analog 
sub-array beamforming without significant reduction in image 
quality. In this way channel number reduction resulting from 
analog sub-array processing is combined with sampling rate 
reduction obtained by FDBF paving the way to real-time low- 
cost 3D imaging system. 

VI. Conclusions 

In this work a solution to one of the major bottlenecks 
in 3D imaging, the amount of sampled data, is introduced. 
The number of samples taken at each transducer element is 
reduced by applying the low-rate sampling scheme presented 
in nil to the individual signals detected by the 2D grid. To 
benefit from the achieved data rate reduction we prove that the 
subsequent processing, namely, 3D beamforming, can be per¬ 
formed directly in frequency. The translation of beamforming 
to the frequency domain allows bypassing oversampling and 
to obtain 4 — 10 fold rate reduction without any assumptions 
on the signal’s structure. 

When signal’s structure is exploited further rate reduction is 
obtained. We prove that 3D beamformed signal obeys an FRI 
model, which allows us to sample and process the signals at 
sub-Nyquist rate while retaining sufficient image quality. 

The performance of the proposed method is verified in terms 
of both LPSF and APSF. It is shown that in accordance with 
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(a) (b) 

Fig. 6: 3D imaging of a phantom of a heart ventricle, (a) displays the time-domain reconstruction of the frame, while (b) 
displays the frequency-domain reconstructed frame, using 1C = B{2 Fourier coefficients of the beamformed signal with 12 
fold rate reduction. 


our expectation it has no effect on LPSF, while the APSF is 
virtually the same when the entire set of Fourier coefficients 
of the beam within its effective bandwidth is computed. For 
sub-Nyquist processing APSF is slightly reduced, however 
when half the beam’s bandwidth is used, the degradation is 
negligible. We also demonstrate the advantage of 3D FDBF 
in the presence of noise. The simulations with noise show 
that low-rate 3D FDBF outperforms not only the time-domain 
processing with the partial grid of elements participating in 
reception, but also the optimal time-domain processing with a 
full grid. 

Finally we incorporate the proposed framework to a com¬ 
mercial imaging system and combine it with analog sub-array 
beamforming, required to adjust the number of elements to the 
number of electronic channels. The results verify that no fur¬ 
ther image degradation is introduced and that our approach can 
be used in conjunction with spatial sub-sampling techniques. 

Our results pave the way for low-cost real-time capability 
crucial for further development of 3D ultrasound imaging. 

Appendix A 

FRI Model of the Beamformed Signal 

We assumed in equation (O that the individual signals obey 
the FRI model. We wish to prove that the beamformed signal 
approximately obeys the FRI model, so that (fTTl) holds. 

In order to show this, we rely on three reasonable assump¬ 
tions. First, we assume that 2(17^1 +| 7 nI+ | 7 ^ ^1) < t/. Such a 
constraint may be forced by applying time-dependent apodiza- 
tion, in such a way that ip{t;0x,0y) is omitted from the delay 
and sum process in (O for t > 2(|7m| + |7n| + l7m,nl)- Sccond, 
we assume that the pulse h{t), transmitted to the medium from 
each of the individual transducer elements and reflected back 
from scatterers in the medium, is compactly supported on the 
interval [0,A). Finally, we assume Act/. Again, such a 
constraint may be forced by applying apodization. 

Using (O, the individual distorted signals in ® are of the 


form 

L 

^x-) ^y) — ^ ^ (^5 ^x-) ^y) (^3) 

1=1 

The resulting signal comprises L pulses, which are distorted 
versions of the pulse h{t). 

Suppose that some of these pulses originated in reflectors 
located off the central beam axis. When we combine the 
individual signals in (| 6 ]) to calculate the beamformed signal, 
these pules will be attenuated due to destructive interfer¬ 
ence. Therefore, when considering the beamformed signal 
^{t;0xj0y), we are concerned only with the pulses which 
originated in reflectors located along the central beam axis. 
For convenience, we assume that all pulses in $23\i satisfy this 
property - those that do not will vanish in ^{t;0x,0y). 

We note that the time of arrival at the {m^n) element, 
is related to the time of arrival at the (mo, no) element 
according to the alignment introduced in (|5]). Thus, we can 
express the delays of the individual signals, in 

terms of t/, as 

tl,m,n ~ 'Tm^nij^V'i^X'i^y)' (24) 

Therefore, we may rewrite (1231) as 

L 

^X") ^y^ — ^ ^ ^xi ^y^i (25) 

1=1 

where we defined h{t;0x,0y) = h{Tm,n{t]^x:^y) — 

Applying our second assumption, we find that the support 
of h{t;0x^0y) is defined by the requirement 

0 ^ ^x-) ^ 7 /) ^xi A. (26) 

It can be shown that the inequalities in (1^ are satisfied for 
t e [t/, t/ + A'), where 

A' = 2 A---, (27) 

tifi + 2A + - 2{-imXe + InVe + 
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(d) (e) (f) 


Fig. 7: Cross-sections of the 3D imaging of a phantom of a heart ventricle, (a)-(c) display the Oy = 0° cross-section, while 
(d)-(f) display the Ox — 0° cross-section. (a),(d) display images acquired using time-domain beamforming, (b), (e) and (c), (f) 
display images acquired with frequency-domain beamforming with 6 and 12 fold rate reduction respectively. 


and we define 

tj,e = \lt‘i+ 4|7m,np - + 7n2/0 + (28) 

Using our first assumption that 2(|7m| + | 7 „| + |7m,nl) < 
and the fact that < 1, we have: 

ti - 2{-fmXe + ^nVe + 7m, 

>tl- 2{\-fmXe\ + hnVel + |7m,n^0l) 
>ti- 2{\^rn\ + |7n| + |7m,nl) 

> 0 . 


Thus, we get A' < 2A, and therefore h{t;0x,0y) = 0 for 
t ^ [7, ^z + 2A). Let us write any t G [7,7+ 2A) as t = 7+7, 
with 0 < 7 < 2A. Then 

7(7 ^x-) Oy^ — + 7? Oy^ Xqji^Yi(ti^ Ox^ 7^)). (29) 

Using our second assumption that A <C t/ and r] < 2A, we 
have T] ^ ti. We then approximate the argument of 7(') in 
(l29l) to first order: 

Xm,n(tl + 7 ? Xj^i^ni^h 0xi (^ 0 ) 

To find the support explicitly, we expand the above inequal- 
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ity. For the left-hand side, we find that 

^y) ^y) — ^l,m,n{^X') ^y) 


(31) 


where 


^l,m,n{^X’! ^y^ — 

1 _ ti - 2{jmXe + InVe + 7m, _ \ 

^ V - HlmXe + JnVe + l^,nZ 0 )tl + M'lrn,n\‘^ J 

(32) 

We now extend our assumption that 2(|7m| + |7n| + 
l7m,nl) ^ and assume that |7m| + |7n| + l7m.nl ^ Hence, 

\ 7 m,n\ = ^J\ 7 m? + | 7 np + | 7 m,nP < l 7 m| + | 7 n| + l 7 m,nl 
ti. Using this assumption, we get cFi^m,n{0x^0y) 1- Replac¬ 
ing T] = t — ti and substituting back to ([29b , results in 

Ox, Oy) ^ /i(t - ti),t G [ti,ti + 2A). 
Substituting back to $25\i and using the fact that h{t — ti) = 0 
for t ^ 2A), we get: 

L 

Oxt Oy) ^ ^ t/). (33) 

1=1 

Finally, plugging this back into (O, 

Ox, Oy) = — ^ ^ Ox, Oy) 

(m,n)eM 

1 ^ 

(m,n)GAt /=1 

i 1 

= (ll,m,nh{t — tl) 

1=1 (m,n)EA4 

L 

= '^bih{t - tl). (34) 

i=i 

Thus, we have shown that the beamformed signal obeys the 
FRI model. 


for t > | 7 m,n|- Assuming the pulse-shape to have a negligible 
support with respect to the penetration depth, A <C T, and 
using the fact that (O holds for all (m, n) G AT, we get, for 
all 1 < / < L: 

min r7„(r;6'^,6»j^). (38) 

{m,n)EA4 

Since denote the arrival times of the echoes to the 

reference element, we can set the upper bound TB{0x,0y) on 

the beamformed signal as: 

TB{0x,0y) = mm Ox, Oy). (39) 

{m,n)EA4 

We are now left to show that TB{0x,0y) < T. This holds 
since we can always find an element (mi,ni) G AT such 
that ^rrn and 7 ^ have opposite signs to that of xq and ye, 
respectively. Furthermore, we note that we can always place 
the reference element (mo, no) such that 7 ^^ m — ^ ^ 

specific choice of (mi,ni) G AT. Thus: 

TB{eM<r-],n,{T-,eM 

T + + \lniye\ 

T 

< T. (40) 

By applying rm,n{t',0x,0y) on both sides of ([39b , we also 
have: 

rm,n{TB{0x,0yy,0x,0y)<T. (41) 


Appendix B 

Beamformed Signal Support 


We consider the FRI model for the individual signals in 
(O. According to our second assumption in Appendix [^ 
h{t) is the known pulse-shape with a support of [0, A) for 
some known A satisfying A <^T. 

We neglect all refiections that reach the (m, n) transducer 
at times greater than T, considering them as noise. Therefore, 
for all 1 < I < L and (m, n) E Ai: 

tl,m,n + A < T. (35) 

Using (El, (E3 and the fact that rm,n{t',0x,0y) is non¬ 
decreasing for t > 0 , we get: 


ti<T-%{T-A;0,,ey), (36) 

with T^]n{'^]0x,0y) being the inverse of Tm,n{t-,0x,(^y) with 
respect to t: 


|7m,ri 


^-1 (-i-.n y \ ^ I (^>^1 _ 

t-{jmXo+lny0+7^,nZo)^ 


( 37 ) 
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